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ABSTRACT 



We report evidence that the asymptotic low-energy power law slope a (below the spectral 
break) of BATSE gamma-ray burst photon spectra evolves with time rather than remaining 
constant. We find that a high degree of positive correlation exists between the time-resolved 
spectral break energy J5 p k and a. In samples of 18 "hard-to-soft" and 12 "tracking" pulses, 
evolution of a was found to correlate with that of the spectral break energy E'pk at the 99.7% 
and 98% confidence levels respectively. We also find that in the flux rise phase of "hard-to-soft'' 
pulses, the moan value of a is often positive and in some bursts the maximum value of a is 
consistent with a value > +1. BATSE burst 3B 910927, for example, has an a max equal to 

■ 1.6 ± 0.3. These findings challenge GRB spectral models in which a must be negative or 
II . remain constant. 
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1. Introduction 

Studies of gamma-ray burst (GRB) spectral evolution have recently begun to uncover trends 

■ which may constrain the emission mechanisms. Earlier reports on spectral evolution focused 



solely on the "hardness" of bursts, measured either by the ratio between two detector channels 



or with more physical variables such as the spectral break or peak power energy E p k (Ford et al. 



1995p which is the maximum of vF v , where v is photon energy and F v is the specific energy flux. 



Such hardness parameters were typically found to either follow a "hard-to-soft" trend (Norris ei 



al. 1986 ), decreasing monotonically while the flux rises and falls, or to "track" the flux during 



GRB pulses ( Golenetskii et al. 1985 ) 



The recent discovery that i? p k often decays exponentially in bright, long, smooth BATSE 
GRB pulses as a function of photon fluence $ provides a new constraint for emission models 



(Liang & Kargatis 1996), and the fact that the decay constant $ often remains fixed from pulse 
to pulse within a single burst hints at a regenerative source rather than a single catastrophic 
event (Meszaros fc Rees 1993). However, that study concentrated only on the evolution of E^. 



To further explore the origin of the spectral break, we begin the analysis of two additional 
parameters in the spectral evolution, the asymptotic low-energy power law slope a below E p ^ 
and the high-energy power law slope [3 above E v y as they are defined in the Band et al. (1993) 
GRB spectral function 
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(a - /3)#o > E, 
(a-[3)E <E, (1) 



where A is the amplitude (in photons sec - 1 cm _2 kev _ ) and Eq = E p \./(2 + a). We note that 
a is not the maximum low energy slope of the GRB function within the detector range, but 
is the asymptotic limit of the slope if extrapolated to arbitrarily low energies. The observed 
values and variability of all three parameters are crucial in evaluating the wide field of proposed 
models of gamma-ray burst emission. For example, many models of the spectral break require a 
to stay constant (e.g. self-absorption) or to have negative values (e.g. — I, Katz 1994 
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T99§), even when E v ^ evolves. In our study, we find that there are hints that /3 decreases over 



the course of some bursts, as suggested by COMPTEL results (Hanlon et al. 1995) and stays 
constant in others. In this letter, though, we focus on the evolution of a and save the discussion 



of (3 for future work (Preece et al. 1997) 



2. Spectral Evolution Patterns 

To determine the evolution of the spectral shape of GRBs, we examine High Energy 
Resolution data collected from the BATSE Large- Area Detectors (LADs) on board the Compton 
Gamma- Ray Observatory ( Fishman et al. 1989] ). We select bursts that have either a BATSE 
3B catalog (Meegan et al. 1996) fluence > 2 x 10 -5 erg cm -2 , peak photon fluxes ~ 10 photons 
s _1 cm~ 2 on the 1024 ms timescale, or a signal-to-noise ratio (SNR) > 7.5 between 25 and 35 
keV with the BATSE Spectral Detector (SD). The inclusion of a few other bursts results in a 
set of 79 bursts. The counts from the LAD most normal to the line of sight of each burst are 
background-subtracted and then binned into time intervals each with a SNR of ~ 45 within 
the 28 keV to 1800 keV range. Employing a non-linear ^-minimization algorithm (Bcvingtor 
~L96|), we fit the Band et al. GRB function to each interval and thus obtain the time evolution 



of the three Band et al. parameters which define the spectral shape. 

Figures 1 through 3 show sample BATSE b ursts (3B) 910807, 9109 27, 911031, 920525, and 
931126, displayed for illustrative purposes. (See Liang fc Kargatis 1996] for E'pk-fluence diagrams 
for 911031, 920525, and 931126.) The spectra show that for these bursts, a generally rises and 
falls with the instantaneous -E p k, though exact correlation between the two parameters is not 
evident. In Figures 1 and 2, /3 stays relatively constant throughout most of the primary pulse, 
while E p k and a both steadily decrease. 

To determine if a does indeed evolve in time in a majority of bursts, we fit a zeroth (M=0) 
and first-order (M=l) polynomial to the a evolution in each burst. Assuming a null hypothesis 
in which a is constant during a burst and the time-resolved values of a are normally distributed 
about the mean, we expect the value A\ 2 = Xm=o — Xm=i t° be distributed as x 2 with 1 degree 
of freedom ( Eadie et al. 1971 ). We calculate for each burst the probability Q of randomly 
drawing a value greater than or equal to A% 2 . We observe that 67 of the 79 bursts have a Q < 
0.05, (which gives a D=0.8 in a K-S test) and 46 have a Q below our acceptable cutoff of 0.001. 
We conclude that a majority of bursts in our sample show evidence for at least a first-order 
trend in a. 

The five sample bursts above suggest that evolution of a mimics that of E p ^. To see if this 
occurs in other bursts, we attempt to disprove the null hypothesis that a is uncorrelated with 
E p k. To test the degree of correlation between a and -E p k, in each of our 79 bursts we compute 
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the Spearman rank correlation r s ( Press et al. 1992] ). For each burst with a positive r s , we 
find the probability P+ of randomly drawing a value of r s that high or higher assuming no 
correlation exists. For each burst with a negative r s , we find the probability P_ of randomly 
drawing a value of r s that low or lower assuming no anti-correlation. The divisions of the bursts 
in this way precludes the inclusion of systematic anti-correlations, which could occur given the 
negative covariance between a and E v \ and the observed shape of our \ 2 minimum contours. 
We next calculate the Kolmogorov-Smirnov (K-S) D statistic between the measured distribution 
of P+ or P- and the distribution one would expect if no correlation or anti-correlation existed. 
We find D=0.45 for the 47 positively correlated bursts. The likelihood of this value, assuming 
no intrinsic correlation, is 2 x 10~ 8 . The bursts showing negative correlation, which suffer 
from systematics described above, are still consistent (likelihood = 0.04) with a non-correlation 
hypothesis. From this we conclude that a positive correlation exists between a and E p \, in at 
least some subset of bursts. 

To determine if this relation exists in "hard-to-soft" or "tracking" pulses, we select 18 pulses 
which we determine to be clearly "hard-to-soft" and 12 pulses which are clearly "tracking" from 
the > 240 pulses within our 79 bursts. Pulses are included in the "hard-to-soft" category if the 
maximum occurs before the flux peak and is greater than E v k at the flux peak by at least 
(TEpk ■ Pulses are "tracking" if the rise and fall of E^ coincides with those of flux to within 1 
time bin (typically ~ \ sec) and if the rise lasts at least 3 time bins. We do not pretend that 
all pulses fall into one of these two categories, but instead treat them as extreme examples in 
a continuum of evolutionary patterns. Following the same analysis described above on these 
smaller populations, we find D=0.46 for "hard-to-soft" pulses and D=0.45 for "tracking" pulses, 
in cases with positive r s . The likelihood of these observed values of r s assuming no intrinsic 
correlations is 0.003 for the "hard-to-soft" pulses and 0.02 for the "tracking" pulses. In contrast, 
while 4 of the 18 "hard-to-soft" pulses and 6 of the 12 "tracking" pulses were anti-correlated, 
the likelihood of these randomly occurring was 0.78 for the "hard-to-soft" cases and 0.29 for the 
"tracking" cases, values consistent with the null hypothesis of no anti-correlation. In Figure 4, 
we compare the cumulative distributions of the 14 "hard-to-soft" and 6 "tracking" pulses which 
are positively correlated to that of the 47 positively correlated bursts. We find that both 
distributions of pulses are similar to the distribution of the bursts which implies an Ep^ — a 
correlation. We conclude from this statistical evidence that for "hard-to-soft" and, with less 
confidence, for "tracking" pulses the asymptotic low-energy power-law slope a evolves in a 
manner similar to -E p k- 



3. Low-Energy Power Index in the Rise Phase of Pulses 

Assuming that a mimics -Epk, it follows that a decreases monotonically for "hard-to-soft" 
pulses, whereas it increases during the rise phase of "tracking" pulses. We compare the averaged 
values of a during the rise phase for these two groups and find that those in "hard-to-soft" 
pulses are significantly higher. While none of the 12 "tracking" pulses has an average a r i sc > 0, 
7 of the 18 "hard-to-soft" pulses had an average a r ; SG > (see Figure 5). A K-S test between 
the two distributions gives a value of D=0.56, implying a probability of 0.014 that these two 
samples were randomly taken from the same distribution. 

We next examine the highest value of a max that occurs in our time-resolved spectra. This 
value serves as a valuable test for GRB emission models. In Figure 6, we provide the distribution 
of otmax found in each of our 79 bursts. Only a few bursts examined so far suggest that their 
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maximum a may be > +1. As indicated in Figure 6, all of the bursts with a max > +1 have 
large statistical uncertainties. The nearly linear decrease of a with respect to time in 3B 910927 
suggests that its relatively high a max of 1.6 ± 0.3, found using data from the LAD most normal 
to the burst, is not merely a statistical fluctuation (see Figure 1). Further examination reveals, 
however, that this burst is still consistent with a < +1 for its duration. In addition, jointly 
fitting the data from the two LADs most normal to the burst reduces a max to 1.03 ± 0.15. 

We also note that fitting 3B 910927 with a broken power law instead of the Band et al. 
GRB function, gives the same linear decrease of the low-energy power law slope 71 with respect 
to time. The fit for the broken power law also has reduced-^ 2 values comparable to those of the 
Band et al. GRB function fit. However, 71 < throughout the burst, a value of the low-energy 
slope lower than that found using the Band et al. GRB function. If the GRB function better 
represents the underlying physics than this difference would be expected. The parameter 71 
measures the effective average slope below E p ^ whereas a measures asymptotic value, allowing 
for the curvature of the exponential function. 



4. Summary and Discussion 

We establish that the asymptotic low-energy power law slope, represented by the Band et al. 
parameter a, evolves with time rather than remaining fixed to its time-integrated value in 58% 
of the bursts in our sample. We find strong evidence that a correlation between the parameters 
.Epk and a exists in the time-resolved spectra of some BATSE gamma-ray bursts and with 
slightly less confidence, we determine that this correlation exists in both "hard-to-soft" and 
"tracking" pulses. We also find that in ~ 40% of the "hard-to-soft" pulses, the average value of 
a during the flux rise phase is > 0, while for "tracking" pulses the average a is always < 0. For 
3B 910927, using data from only the LAD receiving the most counts, we determine a maximum 
value ofa = 1.6±0.3. However, we cannot yet prove that a max > 1 in any burst examined so 
far due to broadness of the % 2 minimum. 

GRB spectral breaks can in principle be caused by synchrotron emission with a low-energy 



cutoff or self-absorption. However, in the former case, a is always < — 1 (Katz 1994, Tavani 



1996| ) with no evolution. Such a low and constant a is inconsistent with many observed 



BATSE bursts. For instance, in 3B 910927, fitting the time bin in which a is maximum with 
an a fixed to — | results in a Q of 1.5 x 10 -11 , much lower the Q=0.35 obtained when a is 
a free parameter. In the case of self-absorption, a could go as high as +1 (thermal) or +1.5 
(nonthermal, power- law) ( |Rybicki fc Lightman 19791) . But again, in such models, a cannot 
evolve with time, only (which would be interpreted as the self-absorption frequency) can. 
Hence, these conventional interpretations of the spectral break of GRB continua can be ruled 



out by our results. Implications of our results on various cosmological scenarios (e.g. Shaviv 



Dar 1996) remain to be be investigated. 



The spectral breaks can also be caused by multiple Compton scattering (Liang fc Kargatii 



1996|) . In this case, the decay of a in "hard-to-soft" pulses can be interpreted as the Thomson 



thinning of a Comptonizing plasma (Liang et al. 1997) and the initial a can in principle go as 
high as +2, because in the limit tt (Thomson depth) — » 00 one would expect a Wien peak. 
However, several factors make it difficult to clearly measure an early low-energy power law 
~ +2 even if the spectral break is related to a Wien peak. The most obvious problem is that 
the highest tt would occur earliest in a "hard-to-soft" pulse, when the flux is the lowest, so 
that fitting a precise spectral model becomes difficult. Another problem is that even if the true 
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GRB spectral break is Wien-like, if one had used a function other than the Band et al. function 
(eg. broken power law) or simply measured the slope within the BATSE range, one could get a 
slope flatter than +2. This is evident in Figure 6, in which the maximum slope for the same set 
of bursts appears to only approach +1 while a max appears to approach +2. This is because the 
exponential curvature depresses the apparent slope relative to the asymptotic power law of the 
Wien function. Also important is that the Band et al. GRB function does not take into account 



the soft X-ray upturn expected from saturated Comptonization of soft photons (Rybicki 



Lightman 197£, Liang et al. 1997). If the lower boundary of the fitting energy window is below 



the relative minimum in the saturated Comptonization photon spectrum (Pozdnyakov, Sobol, 



fc Syunyaev 1983| ), any fitted, low-energy power law, such as the Band et al. GRB function, 
will be flatter than the true slope for the Wien peak. Preliminary results show that moving the 
lower energy cutoff of the fitting region allows one to get a higher a. However, the uncertainty 
in a increases when reducing the size of the fitting window and thus the higher value of a may 
be misleading. Evidence for the X-ray upturns in the low-energy spectra have been found by 
Preece et al. (1996) who found positive residuals between the BATSE data and their fitted Band 
et al. GRB functions in many bursts. However, further analysis of time-resolved, low-energy 
GRB spectra is needed before a model involving saturated Comptonization can be tested. 
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Fig. 1. — Evolution of the Band et al. GRB spectral function for 3B 910927. Each line is marked 
with the time (in s) corresponding to the beginning of the time bin. Note that a typical statistical 
/? error 073 « 0.4. The nearly linear decay of a throughout this burst suggests that the early high 
values of a arc not statistical fluctuations. However, this burst is still consistent with a < +1. 
[Upper inset] Evolution of Ep-^ (squares, logarithmic scale) and photon flux (histogram, linear 
scale) with respect to fluence. [Lower inset] Evolution of a (circles) and photon flux (histogram) 
with respect to time. Error bars represent la confidence level. 

Fig. 2. — Evolution of the Band et al. GRB spectral function of the first pulse in 3B 910807. 
Note that the typical statistical 073 ~ 0.4. See Figure 1 caption for description of plots. 

Fig. 3. — Band et al. a (circles), linear E p k (diamonds) and photon flux (histogram) evolution in 
time for 3B 911031 [GRB 973], 3B 920525 [GRB 1625], and 3B 931126 [GRB 2661]. These three 
bursts suggest that a evolves in a manner similar to E p -^. Error bars represent lcr confidence 
region. The vertical dimensions of diamonds represents la confidence region for and the 
horizontal dimensions represent the durations of the time bins. (See Liang & Kargatis 1996 for 
i? p k vs. fluence of these bursts with a and (3 fixed to the time-integrated value.) 

Fig. 4. — Cumulative distribution of the probability P of randomly drawing a value > t assuming 
no correlation between E^ and a for 14 "hard-to-soft" pulses (squares), 6 "tracking" pulses 
(diamonds), and 47 bursts (circles) with positive correlation. Applying a K-S test, we find 
the probabilities of getting these distributions are, respectively, 0.003, 0.02, and 2 x 10 -8 . The 
probabilities (from a K-S test) that the positively-correlated "hard-to-soft" and "tracking" pulses 
represent samples from the same population from which the 47 bursts are taken are, respectively, 
0.21 and 0.55. 

Fig. 5. — Mean value of Band et al. parameter a for the rise phase of 18 "hard-to-soft" pulses 
(solid) and 12 "tracking" (dashed) pulses. Note that only "hard-to-soft" pulses have mean values 
of a r i S o > 0. The probability of these two samples originating from the same population is 0.014 
as determined by the K-S test. 

Fig. 6. — Cumulative distribution of the maximum value of the Band et al. parameter a (circles) 
and the maximum spectral slope (squares) within the fitted energy region for 79 BATSE bursts. 
The slope (in log-log space) of the Band et al. GRB function equals a — E/E when E < {a—(3)E 
and equals (3 when E > (a—/3)E , where E = E p k/(2+a). Error bars represent the la confidence 
region. Note that slope max represents the maximum slope measured in each burst and does not 
necessarily occur at the time when the Band et al. function a is maximum. Also note that these 
are two separate distributions and the ith a max is not necessarily from the same burst as the zth 
slope max . 
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